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ABSTRACT 



The absorption coefficient for canpressional waves in a water- 
saturated kaolinite clay was ineasured in the frequency range 100 to 
150 kHz xasing a pulse technique with spherically diverging waves. The 
measurement also permitted determination of the sound speed, V^. For 
one condition of the clay with a 77% porosi^^f the sound speed was 
1422 m/sec. The absorption coefficient, a, was found to be frequency 
dependent, satisfying the empirical relationship a = 0.033 f dB/m, f 
being the frequency in kHz. No dispersion of wave speed with frequency 
was observed. These results were used together with the results of 
dynamic shear (rigidity) modulus measurements made with a torsional wave 
viscoelastcmeter to calculate the corrplex Lame constants of a viscoelastic 
model. 
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I. INTRDDUCriON 



For seme problems of sound propagation in the ocean, the ocean bottem 
can be satisfactorily described using the Rayleigh model. Mackenzie [1] 
showed that the Rayleigh model can be applied to cases vhere sound was 
obliquely reflected at grazing angles (measured between the ray and the 
plane of the bottom) greater than the critical angle. However, at 
smaller grazing angles the mcdel can not accoxant adequately for bottom 
losses. In those cases where the attenuation of sound in the sediment 
becones inportant, a more realistic viscoelastic model shoxild be 
employed vhich considers the effects of rigidity and absorption. 

The carplex Laite constants (X/ y, X' and y') of the viscoelastic 
model are theoretically related to four observable quantities, the speeds 
of propagation of cotpressional and shear waves and the attenuation of 
these waves. Ihe effectiveness of the viscoelastic model is therefore 
governed by the precision and accuracy of the values obtained for the 
ootiplex viscoelastic moduli. The moduli can be calculated through measxrre- 
ments of the acoustic properties mentioned above. Measurement of shear 
wave propagation by direct methods is the most difficult of the observations 
to make, particularly in soft sediments, because of its high attenuation 
and the physical prcblom of generating transverse wave motion. One pro- 
cedure to obtain indirectly the shear wave speed is the measurement of 
the speed of a surface wave (Stoneley wave or Rayleigh wave) pjxjpagating 
along the sediment-water interface. When the speed and attenuation of 
the Stoneli^wave is determined, shear wave speed and attenuation constants 
of the sediment can be detennined [2] . Other prerequisites for this lype 



7 



of calculation require the knowledge of corpressional wave speed and 
attenuation and the density of the sediment. Hamilton et al. [3] 
provide a sarrple calculation illustrating this procedure. 

Numerous researchers have conducted tests to determine the corrpres- 
sional wave speed and attenuation ard, to a lesser extent, surface wave 
speed and attenuation of sedirrents either in situ or in the laboratory. 
Urick. [4, 5] .developed expressions for ccnpressional wave speed and 
absorption in kaolinite suspensions. Urick and Jtoent [6] also outlined 
a theory for a ccitplex propagation constant viiose real part yielded a 
wave speed and whose imaginary part yielded an absorption coefficient. 

Laughton [7] conducted laboratory experiments on ocean sediments, 
measuring coipressional and shear wTave speeds only. He was able to make 
direct shear measurements on coirpacted sediments, but not on unconpacted 
material because of its low rigidity. 

Hampton [8] reported on the results of measuranents on the acoustic 
properties of sediments, including conpressional wave speed and attenuation 
measurements of dispersed and landispersed kaolinite sediments. 

Hamilton and his associates at the Naval Undersea Research and 
Developrent Center have conducted both Stoneley and cotpressional wave 
measurements. Their efforts culminated with in situ measurements of both 
waves from the research s\±mersible DEEPSTAR 4000 at stations close 
together and in the same sedimentary environment. With the data, Hamilton 
was able to calculate both elastic and viscoelastic constants for the 
sediments. They were prevented frcm making both sets of measurements 
during the same dive by the quantity of unique equipment required for each 
measurement procedure. 
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Ihis problem may be reduced use of a torsional wave visco- 
elastoneter developed at the Naval Postgraduate School ty Hutchins [9] , 
Cohen [10] and Bieda [11] . The instrument measxores the loading effects 
caused by the rigidity of a sediment on a torsionally vibrating rod. 

These effects can be analyzed to calculate the real and imaginary parts 
of the shear (rigidity) modulus. Determination of the shear nodulus by 
this method could possibly provide a better method to determine shear 
wave properties than by conducting Stoneley wave speed and attenuation 
measurements. 

This experiment was part of a process for validation of the visco- 
elastcmeter measurement method. Since the viscoelastic parameters 
also relate to the soxnd absorption, it is necessary to carefully 
measure the sound speed and absorption in sediment. Octtparison of the 
results of these measurements with those obtained by the viscoelastcmeter 
can give an indication of the validity of the latter method. 

A cotpressional sound wave was propagated by an electrostatic source 
throui^ a water-saturated Jcaolinite sediment. The signal received by a 
hydrophone was recorded and analyzed to obtain values of the cortpressional 
wave speed and absorption. Viscoelastcmeter measurements were then made 
on a core sample of the sediment in vhich the sound wave experiment was 
conducted. 



9 



II. TEffiOKiT 



A. (XlMPRESSIONAL WAVES 

For spherically diverging vaves/ the equation for the soxand pressure 
level (SPL) in decibels can be vnritten [12] : 

SPL = SL - 20 lo?2.0 ^ ~ ^ 

vhere SL is the source level, r is the range and a is the absorption 
coefficient. The logarithmic term represents the spherical spreading loss 
in decibels. By rearranging terms, the relation 

SPL + 20 log^O ^ “ SL - ar (lA) 

provides an equation vhich is convenient for plotting experimental data on 
linear graph paper. The slope of the curve of Equation (lA) should be a 
straight line with units of decibels/unit length, the absorption coeffi- 
cient. Ihe intercept of this curve at a unit distance is the value of 
the source level in decibels. Equation (1) assumes a harogeneous and 
isotropic medium and the absence of reflecting surfaces. For a finite 
sized source, the inverse square law divergence applies only at distances 
from the source that are large in conparison with the radius of the piston 
source. Due to interference effects, the axial pressure near the surface 
of a piston undergoes wide variations with changes in range, until the 
distance is greater than d^/2X, v^ere d is the dianeter of the piston and 
X is the wavelength [13] . 
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B. TEE VISCOELASTIC MCOEL 



Pollcwing the developtnent used by Handlton et al. [3] , the viscoelastic 
model uses coiplex Lame constants (y + iy') and (X + iX‘) to account for 
absorption. For siirple harmonic motion/ the time factor e^“ is used. 
Compressional and vertically polarized shear waves traveling in the sedi- 
ment can be described in terms of displacement potentials and ip. For 
a viscoelastic solid, these displacement potentials must satisfy the 
following wave equations: 



P 




4- = (X + 2y)v2<j, + (X' + 2y') |t 



( 2 ) 



P 




Ip = uV^ip + y ' 




(3) 



where p is the density of the solid. Bucker [14] has solved these 
equations in order to determine the oca^plex laire constants. The resulting 
equations are: 



1/V^ - l/t2 = p/(x + 2y) (1 + q^) 


(4) 


2Aptp = pqp/(l + q|) (X + 2y) 


(5) 


1A| - l/t| = p/y(l + q|) 


(6) 




(7) 



where V and V 
P s 

speeds and 



are, respectively, the ccmpressional and shear wave 



X' + 2y' 



Here, t and t 
P s 

shear waves by: 



are related to the absorption of the oonpressional and 



Sp = -20 (log^Q e) Wt^) (9) 

a^ = -20 (logj^Q e) (oi/tg) 



In this experiment, y and y* were detennined froti measurements by the 
visooelastoneter, and the cxarpressional wave speed and absorption were 
detennined by direct measurement. By using the viscoelastic equations, 

X and X' can be calculated. Values for the expected shear wave speed and 
absorption can also be calculated. 

Assumptions and approximations upon virLch the model is based are 
listed below: 

1. the use of cotiplex lame^ constants to account for attenuation does 
not specif the mechanisms involved; 

2. the dilatational or bulk viscosity is set equal to zero, i.e., 

X' + 2y'/3 = 0; 

3. the carplex Lame constants are independent of frequency, implying 
that the real parts, X and y, governing the cotpressional and shear wave 
speeds, and the imaginary parts, X' and y' , governing energy losses are 
also independent of frequency (most observers, with the notable exception 
of Hampton [8] , have not doserved any dispersion with frequency) ; 

4. absorption in decibels/unit length is assxmed to increase linearly 
with frequency, a = kf (Cole 115] and Wood and Weston [16] have found an 
approximate linear dependence in natural sediments) . 
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C. THE VIS(XIEI^TCME:TER 



The description, theory, and operation of the viscoelastometer have 
been thoroughly reported by Lasswell [17] , Bieda [11] , Walsh [18] and 
others. The viscoelastcmeter, as shown in Fig. 1, consists of a rod- 
shaped torsional irechanical oscillator vhich is excited at one of its 
.mechanical resonance modes by a torsional piezo-electric crystal located 
at its center, setting up a standing wave along the axis of the system. 
The torsional motion of the rod surface in contact with another medium 
generates a shear wave viiich propagates radially outward from the rod. 
The loading effects of the medium, such as a soft clay, cause a measure- 
able change in the mechanical damping and resonance frequency of the 
viscoelastcmeter. Electrical measixrements of these changes can be 
related to the conplex rigidity modulus through a procedure vhich uses 
Newtonian fluids for calibration of the device. 
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III. BQUIPMEnOT AND PPOCEDURES 



A. THE SEDIMENT TANK 

An cirtificial ocean sediment was formed by a saturated kaolinite-water 
mixture contained in 4 ft by 8 ft by 3 ft steel tank (Fig. 2) . The 
tank rests on 9 inches of sand and is surrounded by 1 ft of sand on the 
Icwer half of the .tank's sidewalls. The sand was installed to act as a 
damping mechanism for reflection of waves from the walls of the tank and 
floor vibrations. Kaolinite clay is a hydrated aluminum disilicate with a 
grain diameter comparable to the clay size particles found in many ocean 
bottcm sediments. It was chosen for this experiment for the additional 
reason that it is comtiercially available in a fairly pure state. The depth 
of the clay was approximately 70 cm. 

Ihe bottcm of the tank is lined with ^inch diameter copper tubing 
vhich can be used for heating the clay with steam. Prior to conducting 
measurements, the clay was covered with water, thoroughly stirred in an 
attenpt to break up any clunps and to homogenize it, and heated to drive 
off any gas bubbles that may have existed in the clay. Haitpton [8] 
Qtphasized the importance that even small quantities of gas can have on 
acoustic properties. 

B. OCMPRESSIONAL WAVE MEASUREMENTS 

A horizontally-oriented beam of carpressional waves was generated in 
the clay approximately 25 to 30 cm belcw the clay-water interface a 
7.1-cm diameter electrostatic transducer similar to one constructed and 
described by Palatini [19] . The transducer was a piston-type source with 
a 1-mil thick aluminum coated I^lar film (a DuPont polyester film) 
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attached to a slightly rou^ened flat plate. Ihe transducer was polarized 
by a 300-v battery. Since a thin layer of air was entrapped between the 
plate and the ^fylar film, the air pockets acted as a spring against vdiich 
the film moved vhen the applied alternating signal voltage gave rise to 
an alternating electrostatic force. The motion was expected to be 
essentially uniform over the sxirface and the response was shown by Palatini 
to be independent of frequency at high frequencies in the .mass-controlled 
region where kd>>l (d is the diameter, k is the propagation constant) . 

The electrical drive to the transducer consisted of a gated sine 
wave, usually eight cycles, which was provided by an oscillator and a tone- 
burst generator. The output frcan the tone-burst generator was passed 
throu^ an additional gate and clipping circuit developed by Cohen [10] 
to reduce leak-through vhen the gate is closed. The driving tone-burst 
was amplified in a power aitplifier and stepped up with a 6:1 transformer, 
passed through a polarizing circuit driven by a 300-v battery and then to 
the electrostatic transducer. The frequency of the signal was controlled 
at the oscillator and monitored on a frequency counter. The voltage 
output could be controlled at both the oscillator and the pcwer amplifier 
and was kept constant at 53.5 v by monitoring it on a vacuum tube 

voltmeter. 

The receiver was a ^inch by ^inch cylindrical barium titanate hydrophone 
mounted on a rack-and-pinion assembly to provide smooth and easy positioning 
of the transducer. The received signal was anplified (usually 60 dB) and 
passed through 60- and 120-Hz filters to eliminate lew frequency hum. High 
frequency noise was reduced by also passing the signal through a lew pass/ 
high pass filter set at 180 and 85 kHz. The signal was fed into the boxcar 
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integrator viiere it was scanned by a slowly itioving gate and reproduced on 
an X-Y recorder. Both the gate and the received signal were monitored on 
the oscilloscope. A typical record of the direct arrival sound pulse is 
shown in Fig. 3. A block diagram of the system is shown in Fig. 4. This 
signal processing method had the advantage that it provided records v^ch 
have enhanced signal to iK>ise ratios and vMch can be iteasured more 
precisely than can photographs of the cathode ray oscilloscope. 

The boxcar integrator synchronously sairpled an input signal with a 
variable width # variable delay gate# which can be slowly scanned across 
the input signal. That signal passed by the gate was averaged by variable 
time constant integrators/ the output of vhich was the average of the 
number of repetitions of the input signal over the gate width interval. 
Because the average value of noise over a number of repetitions was zero, 
an inprovement in signal to noise ratio occurred. When the gate was 
scanned across the return signal from the transducer, the waveform was 
reproduced on the X-Y recorder. 

Ihe signals were analyzed to obtain both ccmpressional wave speed and 
absorption. Ihe gate on the boxcar integrator was calibrated with a 
time-mark generator at 10-ysec intervals from 0% to 100% of the boxcar's 
scan/delay scales. The fixed delay (16 ysec) between the trigger of the 
transducer signal and the 0% value on the scan/delay scales was determined. 
Ihe tine of the first arrival of the signal was then read off the waveform 
trace and the delay time was added to get the total travel time. This was 
divided into the horizontal range to determine the conpressional wave 
speed. 
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The abscarption coefficient was determined by a cortparison of the 
anplitudes measured in millivolts at various ranges. During the course 
of the experiment/ it became clear that/ in addition to the direct path 
signal/ there was seme interference caused by possible side wall reflec- 
tions/ sediment-^ter interface reflections and other undetermined 
multipath transmissions. For this reason/ only the first two and one- 
quarter cycles were .used for analysis since they appeared to be free of 

any such interference. The airplitudes in millivolts of the first three 
\ 

positive peaks and the first two negative peaks were measured. Then the 
values for four half-cycles were determined by adding the anplitudes of 
the first positive and first negative peaks / the first negative and 
second positive peaks / and so forth. The values obtained for the four 
segments were then averaged/ and the average was used as a representative 
amplitude for the given frequency and range. The level of this voltage 
was calculated in decibels relative to an arbitrary value of 1 mv by the 
equation: 

SPL = 20 log^^Q (voltage/ rav) 

The reference distance for the spherically divergent term was chosen as 
1 an also for the reason of convenience. Equation (lA) was entered with 
the SPL values out to a range of 1 m and plotted as illustrated in 
Fig. 5. The slope theoretically should be a straight line/ the value of 
v^ch is the absorption coefficient in decibels/ireter. Since scatter 
did exist/ the slope was determined by a least squares fit. A linear 
regression program was run on a Hewlett-Packard 9100A calculator to 
determine the best fit by minimizing the sum of the squares of the 
deviations of the data points fran the straight line. 
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C. THE VISCXEaLASTCMETER 

The visooelastoneter used to provide shear modulus data was the same 
type of instrument used by Bieda [11] and Lasswell [17] . It is 10.5 
inches long with a ^inch diameter cylindrical rod vtose torsional 
oscillation was controlled by a hollcw barium titanate transducer mounted 
in the center of the rod. The transducer was 1.5 inches long and ^inch 
in diameter. It was constructed from an axially-polarized cylinder viiich 
was cut into two parts along the longitudinal axis and rejoined after the 
inversion of one half and the insertion of two electrode grids. The 
transducer was wrapped in fiberglass thread under tension and covered with 
an epoxy resin to improve the contact at the joints. The rod was construc- 
ted frcm two sections of a constant modulus alloy, Nl-Span-C. These parts 
were heat treated and highly polished. 

The visooelastoneter was supported by claitping at the center of the 
transducer losing a tubular steel frairework viiich housed the electrical 
leads to the transducer. In this way the clamp did not interfere with 
the functioning of the viscoelastoreter since the transducer center was a 
node of torsional motion for all the resonant modes. The electrical 
connections were kept watertight by the application of a silicone rubber 
coating. 

The rod was first calibrated in Newtonian fluids of known viscosity. 

The lesonant frequency and electrical resistance were then measured in 
air and in the artificial sediment in order to obtain the changes in 
resonant frequency and electrical resistance required to calculate the 
rigidity. A saitple core of the sediment in vdiich the sound wave experi- 
nent was conducted was obtained and placed in a water bath maintained at 
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a terrperature of 20 C. Three sets of measureanents were taken on the 
core approxiraately 12 hr apart vMch were used by Lieutenant J, Morgan 
to calculate the real and imaginary parts of the rigidity. 



D, WET EENSIIY AND POROSITY 

The wet density is the weight per unit volume of the water-saturated 

sediment. A known volxame of the sediment was obtained by the use of an 

instrument fabricated by the Physics Shop at the Naval -Postgraduate School 

and shown in Fig. 6. This sediment sampler had a known volume of 24.68 
3 

C3tn . The sampler could be lowered to the desired depth m the sediment 
with the only disturbance to the sediment being along the walls of the 
cylinder. When in place, top and bottom plates were rotated across the 
cylinder, trapping the sediment enclosed. The bottom plate contained a 
rubber gasket to prevent leakage of water vdiile the sampler was in the 
upright position. The sample was emptied into a weighing dish and 
weighed. The ratio of the wet weight divided by the volvcte equaled the 
wet density. After being heated in an oven at 105 C for 24 hr/ the 
sample was weighed again, with the weight loss assumed to be entirely due 
to evaporation of the entrapped water. Assuming that 1 g of water occupies 
1 on^, the difference in weight in grams between the wet and dry sample 
was equal to the volume of voids in cubic centimeters in the sample. 
Porosity was then determined by the relationship: 



Porosity 



Volume of Voids 
Total Volume 



Hampton [8] used the term volume concentration vAiich is related to porosity 
by: 

Volume Concentration = 1 - Porosity. 
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IV. RESULTS AND DISCUSSION 



Table I shows the results obtained from the cctipressional wave speed 
measurements in the tank. The speed of sound in water was cxanputed 
using Wilson's formula [12] . The horizontal range measurements were 
accurate to within 0.2 on; the determination of arrival time was accurate 

I 

to withiri'S psec. ' The "calculations -of tlie cortpressional wave speed were 
then accurate to within ±1.5%. The ratio of the sound speed in the 
sediment to that in water is sonev^at lower than would be expected. 

Cohen [10] and Walsh [18] both observed ratios greater than 0.98 for 
kaolinite-water sediments with approximately the same porosity. Data 
collected by Hamilton [20] for shallow water sediments off San Diego 
gave a sound speed ratio of 0.97 or greater. 

A possible explanation for this difference was the existence of seme 
gas bubbles in the sediment. Another indication of this possibility was 
the observation of seme acoustic reverberation in the received sound 
pulse following the direct arrival. Since the reverberation was much 
smaller than the direct arrival, the scattering cross-section per unit 
volume of the inhomogeneities vbich caused the reverberation must also 
have been small. 

The absorption coefficient was determined using Equation (lA) and is 
shown in Table I for the frequencies used in the experiment. Fitting 
these points with a linear equation in frequency resulted in an absorption 
coefficient 

a = kf dB/Vn 
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-2 

with k = 3.3 X 10 , and viiere f is the frequency in kHz. This result 

agrees well with studies made by Hamilton [21] . He predicted a value for 

k of between 0.02 and 0.06 for an artificial sediment with porosity in the 

70 to 80% range. Harrpton [8] repo3rted a value for the absorption of 

approximately 4 dB/ft (approximately 13 dB/in) for kaolinite of nearly 

the same porosity at 100 kHz with a frequency dependence of f * . Wood 

_2 

and Weston [16] determined an absorption coefficient equal to 2 x 10 
dB/ft-kHz (approximately 6.6 dB/k at 100 kHz) in mud that showed a linear 
frequency dependence. 

One of the original objectives of this research was to determine both 
the ccmpressional and shear wave speeds and absorption in the sediment. 
Atterrpts to duplicate Walsh's [18] noasurements of surface waves were in 
vain. Although basically the same equipment and procedures were used, a 
surface wave could not be observed in this sediment. Since a direct cati- 
parison of the results obtained from applying the viscoelastic equations 
to those obtained from the viscoelastometer could not new be made, it was 
decided to assume a value for the shear wave speed in order to make a 
rou^ comparison of the twe methods. 

Values for the real and imaginary parts of the shear imodulus 
(li and y') as computed from viscoelastometer measurements are listed in 
Table II for one of its mechanical resonance iriodes. The values are the 
same order of magnitude as those obtained by Cohen [10] in kaolinite. 
Although Walsh [18] had arrived at a value for the shear wave speed (V^) 
of approximately 23 in/sec, both Cohen's values and the measurements 
reported here indicate that a ^ear wave speed of less than 5 m/sec is 
more realistic for the very soft sedin^t being studied. Values for 
y and y' were calculated using the following viscoelastic equations [3] 
and experimental data: 
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, = pv| 

. 8 « 686ck) 

P 












3q_(X + 2y) 



p = 1.386 g/an^ 

V = 5 X 10^ ati/sec (assumed value) 

0 ) = 2tt(1.56 X 10** Hz) where 15.6 kHz 
was one of the resonant modes of 
the viscoelastometer. 



a = 5.15 X 10”^ dB/an (extrapolated 
from 100 kHz to 15.6 kHz) 

V = 1.422 X 10^ an/sec. 

P 

Table II ccsorpares the results of the two methods. The value of p was 
5 2 

equal to 3.5 x 10 dyne/cm from theory, the same order of magnitude as 

that Treasured by the viscoelastcrreter. The calculated value of p' was 
5 2 

360 X 10 dyne/cm , two orders of irragnitude greater than the viscoelasto- 
neter ireasurement. Examination of the paraireters in the viscoelastic 
equations reveals that the determination of p* is relatively insensitive 
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to changes in the absorption coefficient or the shear wave speed \dien 

V is small. It should be noted that the viscoelastometer results used 
s 

here are approximate since they are based on a very limited number of 
cbservations . However/ they are corparable to observations reported by 
Cohen [10] for a very porous kaolinite mixture. 
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V. EXPERIMENTAL PRCBLEMS 



The sediment, even after being thoroughly stirred, lost its hcmo- 
geneous state in a matter of hours. Settling and (xarpaction took place 
v^ch produced distinct vertical gradients. By lowering one's arm through 
the clay, a definite change in the sediment rigidity could be felt, ranging 
fran a nearly fluid sliorry at the surface to a carpactednmss at the 
bottom. The stratification caused by the changes in density and porosity 
with depth was neither spatially nor terrporally uniform. Horizontal 
gradients, viiile not as pronounced as the vertical gradients, also were 
present. The presence of small clumps of clay caused further discon- 
tinuities in the gradients. Pronounced and changing stratification and 
gradients then did not allow an ideal medium for the propagation of 
sound waves and coiprised a major problem area. 

After the sediment had been reasonably outgassed and homogenized, 
initial atterrpts to obtain conpressional wave measurements were thwarted 
by equipnent problems. Among those were the shorting out of the trans- 
ducer by water leaking into it and low signal-to-noise ratios caused by 
high frequency noise probably generated internally by sane of the 
equipnent. By the time these and numerous minor procedural and equip- 
ment problems were corrected, the sediment had tiire to develop stratifi- 
cation and gradients. Attempts to determine sound absorption coefficients 
shewed that the medium no longer exhibited the inverse square law for 
spherical spreading and was therefore unsatisfactory for further experiments. 
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The sediment was again heated and thoroughly stirred. With increased 
porosity and decreased density, the sound speed decreased as expected [22] , 
but once again the value was low judging by the ratio of the sediment to 
water sound speeds. Some revei±>eration due to scattering fran bubbles 
or other inhomogeneities was c±>served. This co\ald have affected the 
observed absorption coefficient, but the effect is believed to have been 
small. 

Failure to obtain a surface wave measurement indicated other problem 
areas. The sediment in this report was the same kaolinite used by 
Walsh [18] . However, these measurements were made with a layer of water 
overlying it while, in Walsh's ej<periments , the sediment was covered 
with a thin film of water and the surficial layers were much stiffer. 

With the ccnbined effects of the water layer and the extensive stirring 
of the sediment shortly before the experiment, the sediment surface was 
extremely soft. Attempts to place geophones on the sediment surface 
resulted in the geophones sinking several centimeters or more below the 
water-sediment interface. Both the soft nature of the sediment surface 
and the weight of the geophones were contributing factors. Attempts to 
keep the geophones on the surface by placing them on a thin plastic 
sheet also failed. The sheet was pushed beneath the surface vhen the 
geophones were placed on it. As a result, no surface waves could be 
observed. 
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VI. CONCLUSICNS AND RECCWMENDATIONS 



An experimental method has been \osed to determine the speed of propa- 
gation and absorption of ccrrpressional waves in a saturated kaolinite- 

water mixture. For one case of the sediment having a porosity of 

3 

approximately 77%, a wet density of 1.386 g/cm and a terrperature of 20 C, 

the. sound speed was determined to be 1422 itv/sec and the absorption coef- 

_2 

ficient described by a = 3.3 x 10 f dB/m vhere f is the frequency in kHz. 

Data taken by the viscoelastometer on a core of the same sediment confirmed 

earlier results reported by Cohen [10] for basically the same kaolinite- 

water mixture. In addition, Bieda [11] and Lasswell [17] took visco- 

elastcmeter measurements in natural marine sedintents with good correlation 

between their results. Therefore, these results seem to indicate that the 

viscoelastcmeter yields reasonable quantitative estimates of the shear 

moduli. Calculation of \i frcm the observed physical properties (V was 

s 

assumed since it could not be deteimaned) also confirms that the visco- 
elastcmeter measurements appear to be reasonably accurate. However, a 
large discrepancy existed when y' was calculated since the viscoelastic 
equations produced a result vhich was two orders of magnitude greater than 
would be expected for this case. The value of y' is strongly influenced 
by the X + 2y term. This in turn is influenced by the difference in and 
Vp. In sediments with small shear wave speeds, no longer is a significant 
percentage of V , and the calculated y' becomes too large. 

ir 

It is reccmrnended that this discrepancy be resolved by further study. 

If the viscoelastcmeter is producing valid results, and previous data 
indicates that it is, then possibly the theory itay be suspect for low values 
of shear wave speed. One asstmption which itay be questioned is that the 
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V 

V' 




bulk viscosity is zero. If other factors are present viiicdi cause energy 
losses in sediments with small shear wave speed, the present theory may 
not totally account for them and thus nay yield the anomalously high 
value for y*. An exanple of ccnputations by Hamilton [3] of viscoelastic 
constants for a sediment with a shear wave speed greater than 100 m/sec 
illustrated that reasonable results for y' can be obtained from sediments 
that exhibit definite rigidity. 

The data collected in this report for the sound speed and absorption 
were inadequate to provide a complete validation of the viscoelastcmeter 
method. It is necessary to ctotain accurate values for both the ccmpres- 
sional and shear wave speeds and the absorption. With those values, the 
ccnrplex Lame constants for the viscoelastic model can be determined. 

When this experimental pi^ocedure can be satisfactorily acccirplished, then 
results of viscoelastometer measurements on the same sediment can be com- 
pared to validate the effectiveness of the viscoelastorreter method. If 
the viscoelastic theory is foxand to be invalid for small shear wave speeds, 
another sediment having more rigidity iray have to be used in the tank, or 
else in situ measurements of natural sediments will have to be taken. If 
the theory is correct, then there is a serious question of the meaning of 
the imaginary part of the modulus determined by the viscoelastcmeter. 

A recent doctoral dissertation by Borcherdt [23] offers a theoretical 
investigation of inhomogeneous body and surface plane waves in a generalized 
viscoelastic half space. Possibly the theoretical processes can be 
examined in this laboratory environment to provide more information on 
viscoelastic propexrties, in particular, the speed and absorption of shear 
waves. 
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(XMPRESSIONAL WAVE AND SEDIMENT PEDPERTIES 
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TABLE II 



COMPARISON OF METHODS FOR DETERMINING SHEAR MODULI 



Frequency 


Theory 


Viscoelastcmeter 




y 


y' 


• y 


y' 


(kHz) 


(10^ (^e/an^) 


(10^ dyne/cm^) 


15.6 
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360 


1.4 


3.9 
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Figure 1. The Viscoelastcmeter 
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Figure 2. Cross-Sectional View of the 
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Transducer-Hydrophone Separation: 50 cm. 

Scan/Delay : 3 0 %/4 5 % . 

16 ysec delay time has been added to time scale. 
Time scale can be linearly interpolated. 
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Figure 4. Block Diagram for Ccitpressional Wave Ifeasurements 
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Figure 5. Plot of SPL + 20 log^Qr vs Pange (r) 
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Figure 6. Sediment Sanpler 
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